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Abstract. - We analyze the noise induced synchronization between a collective variable char- 
acterizing a complex system with a finite number of interacting bistable units and time periodic 
driving forces. A random phase process associated to the collective stochastic variable is defined. 
Its average phase frequency and average phase diffusion are used to characterize the phenomenon. 
Our analysis is based on numerical solutions of the corresponding set of Langevin equations. 



Introduction. — Synchronization phenomena in non- 
linear systems have been widely studied in very different 
contexts [1-3]. In complex nonlinear autonomous systems 
with many interacting subunits, nonlinear dynamics might 
lead to self-synchronization between the individual ele- 
ments. Noisy environments and externally applied time 
periodic forces might also influence synchronization. One 
might be interested on how noise and driving forces af- 
fect the synchronization between the different subunits, 
or the focus of interest might be the noise induced forced 
synchronization to the external driving force of a suitably 
defined single collective variable characterizing the system 
as a whole. 

In the case of systems characterized by just a sin- 
gle noisy dynamical variable subject to the action of a 
deterministic external time-periodic signal, noise-induced 
forced synchronization has been amply studied. With an 
appropriate definition of a random phase process asso- 
ciated to the response of the system, theories of noise- 
induced forced synchronization have been developed and 
analytical expressions for the quantities characterizing the 
forced synchronization mechanism have been obtained 
[4-9]. The theory of forced synchronization in a single 
variable system has also been extended to the quantum 
regime of incoherent tunneling transitions in driven dis- 
sipative systems [10]. Synchronization effects in complex 
systems formed by arrays of excitable systems and stochas- 
tic resonators have been analyzed in [3,11]. 

In this work, we consider a model describing a finite set 
of N interacting bistable subsystems, each of them char- 



acterized by a single degree of freedom Xi (i = 1, . . . , N), 
whose dynamics is governed by the Langevin equations 
[12,13] " 



Xi — Xi 



i)+Si{t)+F(t). (1) 



Here, the £j (t) 's are Gaussian white noises with zero av- 
erage and (£i(i)£j(s)) — 2D6ijS(t — s), 8 is the parameter 
defining the strength of the interaction between subsys- 
tems, and F(t) — F(t + T) is an external driving force 
of period T. The amplitude of the driving force is sup- 
posed to be large so that a linear response approximation 
of the dynamics is not adequate. Nonetheless, we are in- 
terested in noise induced synchronization, so the driving 
amplitude is subthreshold in the sense that, in the absence 
of noise, the external driving by itself can not induce the 
phenomenon. 

We focus our interest on a single global variable, rather 
than on the variable characterizing an individual subsys- 
tem. In a recent work, we have analyzed the phenomenon 
of Stochastic Resonance associated to the global variable 
of the same system [13, 14]. In this work, we analyze the 
phase process associated to the noise induced jumps of the 
global variable between its dynamical attractors. Noise in- 
duced forced synchronization between the global variable 
and the driving force is characterized by the behavior of 
the average phase frequency and the average phase diffu- 
sion constant with the noise strength. 
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The random phase process, its average phase fre- 
quency and the diffusion constant. The variable 
of interest is the collective output S(t) defined as 



S(t) 



1 N 



(2) 



In the limit N — > oo, an asymptotically valid Langevin 
equation for S(t) can be obtained from Eq. (1) [15]. For fi- 
nite systems the validity of this Langevin equation is ques- 
tionable and we rely on the solution of the Langevin equa- 
tion for each xi to obtain information about the global pro- 
cess S(t). For low noise strengths and driving forces with 
sufficiently large periods, our numerical simulations show 
that a random trajectory of S(t) contains essentially small 
fluctuations around two values (attractors) and random, 
sporadic transitions between them. For each realization of 
the noise term, we then introduce a random phase process, 
4>(i), associated to the stochastic variable S(t) as follows. 
We refer to a "jump" of S(t) along a trajectory, when 
a very large fluctuation takes the S(t) trajectory from a 
value near an attractor to a value in the neighborhood 
of the other attractor. We count N { a '(t), the number of 
jumps in the a trajectory of the process S(t) within the 
interval (0, t] . A trajectory of the phase process is then 
constructed as 

(j>^ (t) =^iV (Q) (i), (3) 

so that <f>(t) increases by 2tt after every two consecutive 
jumps. 

The first two moments of the phase process are esti- 
mated as 



M 

w*)> = ^E* (a) (*) 



(4) 



a=l 



i m 

= i<E[* ( "(') 



M 



.a=l 



(5) 



where, M. is the number of generated random trajectories 
(typically, 3000 trajectories for the results presented in 
this work). 

The instantaneous phase frequency is easily determined 
from the time derivative of {(j>{t)). After a sufficiently large 
number of periods of the driving force, n, the system for- 
gets its initial preparation, but the instantaneous phase 
frequency is still a function of time. Then, we define a 
cycle average phase frequency il p h by averaging the in- 
stantaneous phase frequency over a period of the external 
driving, [7, 10] 



1 [W? M d(<j>(t)) (<j> [(„ + - W.nT)) 

nT "' 

(6) 



%h J- [ <H 



Similarly, the cycle average phase diffusion coefficient 
is evaluated from the instantaneous slope of the variance 
v(t) as [7,10], 



^ 1 



D 



P h 



T J nT dt 



v [(n + 1)T] - v (nT) 



(7) 

In previous works, approximate anlytical expressions for 
these two quantities have been derived for the N = 1 
problem in the classical [4,7,9,16] and quantum cases [10]. 
Those expressions can not be applied to the collective vari- 
able of an TV-particle problem, as a closed Langevin or 
Fokker-Planck equation for S(t) for a finite size system 
does not exist. Thus, we will rely on numerical solutions 
of the Langevin equations, Eq. (1), for the estimation of 
the average phase frequency and diffusion coefficients in 
Eqs. (6) and (7). The numerical method used to solve the 
Langevin dynamics has been detailed in [17]. 

Results. — We have analyzed the forced syncroniza- 
tion phenomenon for the collective variable S(t) for sev- 
eral values of the coupling strength, 6, the noise strength, 
D, and two types of periodic forces: sinusoidal forces 
(F(t) = A cos Fit) and rectangular forces (F(t) = A 
(F(t) = -A) if t e [nT/2, (n + l)T/2) with n even (odd)). 

In Fig. 1 we depict the behavior of the average phase fre- 
quency, f2 p h, and diffusion constant, £> p h, obtained from 
numerical simulations. The external driving force is taken 
to be sinusoidal with A — 0.3 and — 0.01. In the ab- 
sence of noise, the variable S(t) does not jump between 
the attractors. Thus, for these parameter values, the ex- 
ternal driving is subthreshold and we have noise induced 
effects. It is clear from the picture that, even for a single 
subsystem, N = 1, forced synchronization exists as there 
is a range of noise values for which f2 p h matches the ex- 
ternal frequency, £1 and the diffusion constant is rather 
small, with a minimum around D w 0.02. As we discused 
previously in [7], these numerical results agree very well 
with the analytical expressions reported in [7] . For a set of 
N = 10 independent subsystems (8 = 0), forced synchro- 
nization for the collective variable is very much enhanced 
with respect to that observed in a single particle variable. 
The matching between the driving frequency and the av- 
erage phase frequency extends over a very large range of 
noise values, for which the corresponding diffusion con- 
stant is very much reduced. It should be pointed out that, 
in this and subsequent figures, the fact that a "range" of 
D values shows up in the plots, for which D p h seems to be 
constant, is due to the lack of machine precision to com- 
pute very small numbers. It is safe to say that for the 
noise values in those ranges, the phase diffusion constant 
is zero, and that the collective variable jumps every half- 
period of the driver without skipping any change in sign 
of the applied force. 

The results for a rectangular input force with an am- 
plitude A = 0.3 and fundamental frequency £1 = 0.01 are 
shown in Fig. 2. The behaviors of £! p h and D p h with the 
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Fig. 1: The behavior of the average phase frequency, fl p h, and 
the phase diffusion constant, D p h, with the noise strength, D, 
for a sinusoidal driving term with A — 0.3 and Q — O.Of and 
several values of the interaction strength 6. 



noise strength for the different sets of parameter values is 
qualitatively similar to the one observed in Fig. 1 for a 
sinusoidal driving. Quantitatively, for a rectangular driv- 
ing, the noise induced forced synchronization is enhanced 
with respect to that observed with a sinusoidal input with 
the same amplitude and fundamental frequency, as indi- 
cated by the wider range of D values for which frequencies 
match. 

We have analyzed the influence of increasing the driving 
fundamental frequency or decreasing the external ampli- 
tude. In Fig. 3 we depict the results obtained for O p h and 
£>ph for a rectangular input with A = 0.1 and £1 = 0.01, 
with N — 10 interacting subsystems with 9 — 0.5. Com- 
parison with the corresponding results in Fig. 1 indicates 
that synchronization for this smaller amplitude driving 
is much weaker than for the larger amplitude case. Not 
only the noise range for which the average phase frequency 
matches the external frequency is reduced, but the dip in 
the phase diffusion constant plot is much less pronounced. 

In the next figure, we plot the results of our simula- 
tions for a rectangular input with A = 0.3 (the same as in 
Fig. 1) and = 0.1, with N = 10 interacting subsystems 
with 9 = 0.5. The increase in the fundamental driving fre- 
quency with respect to the value used in Fig. 1 leads also 
to a reduction of the noise range for which noise induced 
forced synchronization is optimal. Nonetheless, the dip in 
the Dph vs. the noise strength seems to be of the same 
order of magnitude as observed for the lower frequency 
case depicted in Fig. 1. 

The weakening of noise induced forced synchronization 
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Fig. 2: The behavior of the average phase frequency, f2 ph , and 
the phase diffusion constant, D pb , with the noise strength, D, 
for a rectangular driving term with A — 0.3 and fundamen- 
tal frequency Q = 0.01 and several values of the interaction 
strength 6. 
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Fig. 3: The behavior of the average phase frequency, fl p h, and 
the phase diffusion constant, D p h, with the noise strength, D, 
for a sinusoidal driving term with A = 0.1 and Q, — 0.01. 
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Fig. 4: The behavior of the average phase frequency, fi p h, and 
the phase diffusion constant, D p h, with the noise strength, D, 
for a rectangular driving term with A = 0.3 and — 0.1. 



as the driving amplitude and/or the driving period de- 
crease is to be expected. Noise induced forced synchro- 
nization is due to the modification by the applied force of 
the intrinsic noisy dynamics of the set of bistable systems. 
Noise itself makes the collective variable jumps between 
attractors with a certain probability distribution. The ex- 
ternal driving greatly modifies this distribution by forcing 
the jumps to be coherent with the periodicity of the exter- 
nal driving. The external force tilts the potential surface 
such that every half-period, jumps towards one attrac- 
tor arc much favored against the ones towards the other 
one. Thus, stronger amplitudes and longer driving peri- 
ods, provide potential surfaces which remain tilted during 
a time interval long enough for the noise induced jumps of 
the subsystems variables to happen almost simultaneously 
every half-period of the driver. 

In conclusion, we have numerically demonstrated the 
phenomenon of noise induced forced synchronization in 
finite arrays of bistable systems. Interactions between the 
subunits do not destroy this coherent effect. Indeed, the 
collective synchronization is much enhanced with respect 
to the one observed for a single bistable unit. Driving 
forces with short periods and small amplitudes tend to 
weaken drastically the phenomenon. 
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